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Abstract. Wc study the hydrogen atom confined to a spherical box with 
impenetrable walls but, unlike earlier pedagogical articles on the subject, we assume 
that the nucleus also moves. We obtain the ground-state energy approximately by 
means of first-order perturbation theory and by a more accurate variational approach. 
We show that it is greater than the one for the case in which the nucleus is clamped 
at the center of the box. Present approach resembles the well-known treatment of the 
helium atom with clamped nucleus. 



1. Introduction 

Confined quantum-mechanical models have proved to be suitable first approximations 
for estimating the effect of pressure on the spectral lines of atoms and molecules or the 
effect of their neighbours in condense media. Several such models have been proposed 
for pedagogical purposes in introductory courses on quantum mechanics [TM20]. Among 
them we mention the quantum bouncer [2ll3ll6] . the harmonic oscillator [8]I121I13II16P18II20] 
and the hydrogen atom [3|9|,[TTl[T2l[Tll[T5l[T71[l9]. Such models have also been useful for 
the discussion of semiclassical approaches [HIE], the variational method p21[T5],[20] and 
perturbation theory [3j[T8j[20] • Regarding the latter approach we mention that the sum 
over states is impractical for the calculation of corrections of order greater than the first 
one [3tfl8|l20]. It is preferable either to integrate the perturbation equations directly [2T] 
or to make use of the hypervirial and Hellmann-Feynman theorems [21|[22]. 

In the case of the harmonic oscillator, most studies refer to the case of a particle 
that moves in a box under the effect of a potential of the form V(x) = k(x — xq) 2 /2 as 
if it were tied to the point x by means of a spring of force constant k [8, 12l[T6 ][T8l[20] . 
If we assume that the Hooke's force is due to the interaction between two particles then 
we have a model similar to the one discussed by Tanner [13J who showed that it is not 
possible to separate the center of mass and internal degrees of freedom in the usual 
way because of the effect of the boundary conditions. Amore and Fernandez [23] have 
recently discussed this problem in greater detail. 
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The usual model for the confined hydrogen atom suffers from the same limitation: 
the nucleus is considered to be clamped somewhere inside the box [7J[9l[IIl[12l[]3][15] ' ^ 
appears to be most interesting to assume that not only the electron but also the nucleus 
moves inside it. The purpose of this paper is to discuss such a model in the simplest 
possible ways. In Sec. |2] we outline the model and write the Schrodinger equation in 
a dimensionless form. In Sec. |3] we obtain the ground-state energy approximately by 
means of first-order perturbation theory and compare it with the one for the clamped- 
nucleus case. In Sec. H] we resort to a somewhat more elaborate variational function. 
Finally, in Sec. El we summarize and discuss the results. 

2. The model 

The Hamiltonian operator for a nonrelativistic hydrogen-like atom is 
H =f + V 

~ h 2 ^ 2 h 2 ^ 2 

T = V V 

2m e e 2m n n 

Ze 2 

V(r) = - -p— (1) 
47re r 

where m e and m n are the masses of the electron and nucleus located at r e and r n with 
charges — e and Ze, respectively, r = |r|, r = r e — r n , eo is the vacuum permittivity and 
V 2 denotes the Laplacian in the coordinates indicated by the subscript. 

In the case of the free atom we separate the motion of the center of mass from the 
internal one by means of a well-known change of variables and obtain 



H — Hcm + H int 

h 2 

2M 



H cm = ~ ttttVcm^ M = m e + m % 



where V 2 and Vq M are the Laplacians for the variables r and t C m = (rn e r e + m n r n )/M, 
respectively, and m is the reduced mass. Thus, we can factor the energy states of the 
free hydrogen atom as ip(r e , r„) = ipcM{ r CM)ipint{ r ) an d solve the Schrodinger equation 
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for H int in the usual way [2UI25]- The eigenfunctions and eigenvalues of H int provide 
all the physical properties of the isolated atom, such as, for example, the spectral lines, 
selection rules, etc. [24^125]. 

If the atom is confined to a spherical box of radius R with impenetrable walls then 
the states should vanish when either r e = R or r n = R and, consequently, the above 
separation is not possible as discussed by Tanner [13] and Amore and Fernandez [23] 
for the harmonic oscillator. The reason is that the variables Tcm and r are unsuitable 
for the boundary conditions that are naturally given in terms of r e and r n . 

The positions of the electron an nucleus in the box are completely determined by 
six variables. We conveniently choose r e , r n and r (the sides of a triangle) plus three 
angles for the orientation of the triangle in space. The S states (those with zero angular 
momentum) depend only on the three radial variables: ib(r e ,r n ,r). In fact, if we take 
into account that 

r e dib r dib 
r e or e r or 

V ih = _ (3) 

r n dr n r dr 

then we realize that ib(r e ,r n ,r) has zero angular momentum: 

(r e x V e + r n x V„) V = ^ZlAjll^ = (4) 

r or 

In order to simplify the calculation we first make the change of variables r' e = r e / R 
and r' n = r n / R that leads to the dimensionless Hamiltonian operator 



7 12 _ !_V7'2 

/>- " 2 ' e 2 

m e , m P Ze 2 R 



The states of this dimensionless system vanish when either r' e — 1 or = 1. From now 
on we omit the primes on the dimensionless quantities but keep in mind that lengths, 
masses and energies are measured in units of R, m e , and h 2 /(m e R 2 ), respectively. For 
example, 1//3 is the nuclear mass in such units. 
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3. Perturbation theory 
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For simplicity we restrict ourselves to the ground state and a small box radius. If A 
is a small parameter then we can try perturbation theory in terms of the unperturbed 
or reference Hamiltonian = Hd(X — 0). The perturbation is therefore given by the 
interaction between the particles H' d = — 1/r. The unperturbed ground state is 

sin(vrr e ) sin(vrr n ) 

<p{r e , r n ) = 2 (6) 

r e r n 

Therefore, the expectation value of Hd with this function gives us the energy of the 
ground state corrected through first order of perturbation theory. Besides, according to 
the variation principle such approximate energy will be an upper bound to the exact 
one [MI25]. 

The calculation is reminiscent of that for the helium atom under the clamped- 
nucleus approximation and we may therefore profit from well-known results. The 
calculation of the expectation value of the kinetic energy is straightforward, and there 
are various ways of calculating the expectation value of 1/r [2UI25]. Here, we resort 
to the expansion of 1/r in terms of Legendre polynomials that leads to the simple 
integral [24J 



(f(r e ,r r 



\2 



dr P dT„ = 16n 2 



l rr, 



(p{r e ,r n ) 2 r e r 2 n dr n dr e 



UQ JO 
1 rl 



+ / (f(r e ,r n ) r e r n dr n dr e 

'0 Jr e 



(7) 



Since the analytical expression is rather cumbersome we just show the numerical result: 

e(A) = + X) - 1.786073167A 

= 4.934802200(/3 + 1) - 1.786073167A (8) 

Computer algebra systems are nowadays available in the science departments of 
most universities because they are invaluable teaching tools. This problem may be 
useful for motivating the students to resort to such software. 
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We can obtain simple analytical expressions by means of the even simpler trial 
function 

p(r e ,r n ) = 30(l-r e )(l-r ft ) (9) 
that leads to a quite similar result 

e(A) = 5(0 + 1)-^ 

= 5(/3 + 1) - 1.785714285A (10) 

It is interesting to compare the results for this model with those for the hydrogen 
atom with the nucleus clamped at the center of the box. If we calculate the expectation 
value of the dimensionless Hamiltonian operator (notice that we use the same units as 
before) 

Hm = - £ (11) 



with the approximate trial function (p(r) = V30(l — r) we obtain 

e*(A)=5- y (12) 

For comparison we also consider the unperturbed ground state 

. , ^sin(7rr) , . 

tp r = (13) 
r 

that leads to the first-order dimensionless energy 

e H (\) = — - 2.437653392A 

= 4.934802200 - 2.437653392A (14) 

Obviously, these results are valid for sufficiently small values of A. After contrasting 
equations ( |T2l and ( fl4l) with more accurate results [2T] we conclude that present first- 
order estimates are acceptable for A < 1. In principle, we may assume that the accuracy 
of present moving-nucleus results are as accurate as the clamped-nucleus ones for A < 1 
. If this is true, then our results suggest that the energy of the moving-nucleus model is 
larger than the clamped-nucleus one (e(A) > e#(A)), at least for sufficiently small box 
radii. The difference does not come mainly from the kinetic energy of the nucleus that is 



The confined hydrogen atom with a moving nucleus 7 
proportional to /3 ~ 1/1836 but from the electron-nucleus interaction. This conclusion is 
consistent with earlier variational results that show that the smallest energy takes place 
when the nucleus is clamped at the center of the sphere and increases as it approaches 
the wall pj)] . Fig. CD shows the approximate energies given by equations (E]) , ( TTOT) , ( |T2l) 
and f|T4|) for A < 5 as well as accurate numerical energies for the clamped nucleus model 
calculated by a straightforward power-series method |21j . 

The critical value of A defined by e#(A c ) = estimated from the first-order 
perturbation energy (fT4"|) A c ~ 2 is about 9% larger than the actual value A c = 
1.835246330 that one easily obtains by means of the method already mentioned above 
or from perturbation theory of greater order [21]. For the moving-nucleus model our 
approximate expressions ([8]) and ffTOl predict A c ~ 2.8 and we expect that its error is of 
comparable magnitude. 



4. Variational method 



We can improve the results of the preceding section by means of the variational method. 
To this end we propose the trial function 

<p(r e , r n , r) = N(a)(l - r e )(l - r n )e~ ar (15) 

where a is a variational parameter and N(a) an appropriate normalization factor. Since 
it has the correct asymptotic behaviour e~ ar for the free atom, we expect it to yield 
accurate energies for large values of A. More precisely, we expect that it yields the exact 
energy of the ground state of the free hydrogen 

ell 

lim — = ; — ~ — (16) 

A^ooA 2 2(1 + 0) 2 K ' 

By means of a straightforward but tedious calculation we prove that 

1 / d 2 2d r 2 -r 2 +r 2 8 2 



2 \^dr 2 r e dr e ^ r e r dr e dr ^ 
13 ( d 2 2 8 r 2 -r 2 + r 2 d 2 



2 \dr^ r n dr n r n r dr n dr * 

l + Bfd 2 2d\ A 
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It follows from r 2 = r 2 + r 2 — 2r e • r n and —r e r n < r e • r n < r e r n that 
[r" e — T n | < r < r e + r n . Therefore, we can calculate all the integrals that are necessary 
for the variational approach by means of the simple recipe [26] 

/ / f(r e ,r n ,r)dT = / / / f(r e ,r n ,r)r e r n rdr e dr n dr 

J J J Jo Jo Jr e —r n 

rl rr n i>r e +r n 

+ / / X r e,r n ,r)r e r n r dr n dr e dr (18) 

JO Jo Jr n -r e 

The calculation of the variational energy is simple but tedious and yields rather 
cumbersome results. In what follows we simply outline the main steps. First, notice 
that 

e(a, A) = = (T) (a) - X (1/r) (a) (19) 

is a linear function of A. Therefore, instead of calculating a (A) from the minimum 
condition de(a, X)/da = we derive an analytical expression for A (a) and a parametric 
expression for the variational energy: e(a,X(a)). One can easily carry out this 
calculation by means of a computer algebra system. 

Fig. |2] shows approximate values of e/A 2 obtained from the trial functions flH]) and 
( TT5l) . We appreciate that there is a satisfactory agreement for A < 1 as expected but 
it seems that our earlier prediction of the critical value A c was not as accurate as we 
believed. The variational function f lT5|) yields the more accurate value A c = 2.262 that 
differs 21% from our earlier estimation. We also realize that the variational energy 
approaches the exact free-atom energy (|T6|) as A increases. 

In passing, we mention that the variational calculation confirms the earlier 
perturbation result that the energy of the moving-nucleus atom is larger than the one 
with the nucleus clamped at origin. 

5. Conclusions 



Tanner [13J proposed a pedagogical discussion of the effect of the boundary conditions 
on the separability of the degrees of freedom of a confined system. However, he did 
not show any result for the one-dimensional harmonic oscillator that he chose as an 
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illustrative example. Later Amore and Fernandez [23] discussed that model in more 
detail. In this article we extended those arguments to the hydrogen atom and carried 
out simple approximate calculations for the ground state by means of straightforward 
first-order perturbation theory and a more accurate variational approach. It has been 
our purpose to show how to do the calculation using well-known techniques already 
applied to the helium atom with the clamped-nucleus approximation. Our analysis 
shows that it is not possible to separate the Schrodinger equation in the usual way in 
terms of internal and center-of-mass coordinates r and yqmi respectively, because the 
boundary conditions are given in terms of the electron and nucleus coordinates r e , and 
r n , respectively. Our approximate results show that the energy is greater when the 
nucleus moves than when it is clamped at the center of the spherical box and that the 
difference does not come mainly from the kinetic energy of the moving nucleus that 
is considerably smaller that the electronic one. Present results based on perturbation 
theory are limited to a small box radius or strong confinement but the variational ones 
are valid for all box sizes and yield the free-atom energy when R — > oo. Although it 
is relatively easy to carry out perturbation calculations of large order for the clamped- 
nucleus model [21], the treatment of the moving-nucleus case is considerably more 
complicated. The variational method appears to be a better choice. It can be improved 
as in the case of the Helium atom by means of a Hylleraas-like expansion [26] for the 

trial function ip(r e ,r n ,r). 
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Figure 1. Results from equations ((8|) (circles), ([T0|) (solid line), (fT2]) (dashed line), 
(fT^jl (squares) and power series (filled circles). 
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Figure 2. Approximate values of 4^ obtained from the trial functions (0) (dashed 
line) and (IT5|) (solid line). 



